Signature of pseudogap formation in the density of states of underdoped cuprates 
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The resonating valence bond spin liquid model for the underdoped cuprates has as an essen- 
tial element, the emergence of a pseudogap. This new energy scale introduces asymmetry in the 
quasiparticle density of states because it is associated with the antiferromagnetic Brillouin zone. 
By contrast, superconductivity develops on the Fermi surface and this largely restores the particle- 
hole symmetry for energies below the superconducting energy gap scale. In the highly underdoped 
regime, these two scales can be separately identified in the density of states and also partial density 
of states for each fixed angle in the Brillouin zone. From the total density of states, we find that the 
pseudogap energy scale manifests itself differently as a function of doping for positive and negative 
bias. Furthermore, we find evidence from recent scanning tunneling spectroscopy data for asymme- 
try in the positive and negative bias of the extracted A(^) which is in qualitative agreement with 
this model. Likewise, the slope of the linear low energy density of states is nearly constant in the 
underdoped regime while it increases significantly with overdoping in agreement with the data. 

PACS numbers: 74.72.-h,74.20.Mn,74.55.+v 



I. INTRODUCTION 



BCS theory modified to account for d-wave symme- 
try of the superconducting order parameter has provided 
a solid basis for a first understanding of the properties 
of the cuprates around optimum and in the overdoped 
regime. However, the underdoped region of the phase 
diagram provides challenges to such a simple approach. 
These systems present many features that have been con- 
sidered anomalous and are not part of simple BCS the- 
ory. While the subject remains controversial and very 
different ideas have been put forward to understand these 
anomalies, a recently developed model by Yang, Rice and 
Zhang (YRZ)i which is based on a spin liquid resonating 
valence bond approach, has had considerable success in 
this direction. It is very different from the preformed pair 
model^ where a superconducting gap forms at a high tem- 
perature, the pseudogap temperature T* , and the super- 
conductivity appears only at a lower temperature Tc as a 
result of the onset of phase coherence. The YRZ model 
has a central element, a new energy scale, the pseudogap 
Apg, which is responsible for changes in the electronic 
structure of the normal state above Tc- The pseudogap 
is distinct from the superconducting gap Age and the su- 
perconducting state is conceived as forming from this new 
normal pseudogap phase which is quite different from an 
ordinary Fermi liquid state. In the YRZ model, the large 
Fermi surface of Fermi liquid theory (FLT) reconstructs 
into hole and electron pockets as a result of the growth 
in pseudogap at doping levels x, in the underdoped re- 
gion of the cuprate phase diagram, with x < Xc where 
Xc is the doping associated with a quantum critical point 
(QCP). 

The model is related but different from other compet- 



ing order proposals such as D-density wave formation-^ii 
and has the desirable property that, in its final form, 
it remains simple and has been successfully applied to 
the calculation of many superconducting properties'^—. 
Besides a self-energy which accounts for the formation 
of a pseudogap on the antiferromagnetic Brillouin zone 
(AFBZ) boundary, the model has Gutzwiller factors 
modifying the underlying band structure parameters. 
These narrow the bands as correlation effects become 
more important. Also, a Gutzwiller factor accounts for 
the loss of coherence which greatly reduces the weight of 
the remaining quasiparticle peak as the incoherent back- 
ground increases. These elements account for the ap- 
proach to the Mott insulating state which is best under- 
stood near half-filling in a localized picture of electron dy- 
namics. The hopping from one site to another is blocked 
by a large Hubbard U which describes the energy cost for 
double occupancy. The pseudogap and AFBZ then plays 
a role similar to a band gap at the Brillouin zone in or- 
dinary band theory but with essential differences. For 
example, in YRZ, the bands are not filled rigidly with 
decreasing doping but instead undergo profound changes 
as the pseudogap increases. 

Among the properties already calculated and com- 
pared with experiment are Raman^:- , specific heatX, op- 
tical properties*, aspects of angular resolved photoemis- 
sion spectroscopy (ARPES) and scanning tunneling spec- 
troscopy (STS)2 including the checkerboard patterni^, 
Andreev tunneling^ ^ and also the penetration depth'^. 
Each of these properties show behaviors which cannot 
be understood in d-wave BCS nor in it extensions to in- 
clude inelastic scattering^"^"^*^, anisotropy^i"— , or strong 
coupling effects rooted in Eliashberg theory2i"— . Among 
the previously considered anomalous properties that are 
now understood are the two distinct gap scales seen in 
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Raman spectra. The peak decreases with decreas- 
ing doping while the Big scale increases insteadS,iS. The 
normalized jump in the specific heat drops rather precip- 
itously as X decreases towards the bottom of the super- 
conducting dome2. New structures are seen in the optical 
self-energy and the two energy scales found in the partial 
optical sum as a function of energy are understood- . The 
rapid drop in the superfluid density at zero temperature 
with decreasing x while the slope of the low tempera- 
ture linear in T law is relatively only weakly changed is 
also showiii^. Encouraged by these successes, here we 
consider the quasiparticle density of states (DOS). 



In Sec. II, we present the formalism of YRZ^ for the 
electron spectral density in the underdoped regime. It in- 
cludes pseudogap formation below a QCP at a doping of 
X = Xc = 0.2. We describe how this new energy scale 
modifies the electronic structure from the usual large 
Fermi surface of FLT to Luttinger pockets. There are 
two energy branches in the theory , with Ej^ giving 
a hole pocket and E'^ an electron pocket, the latter only 
in a restricted doping range just below Xc- For very un- 
derdoped samples, only the hole pocket remains. In Sec. 
Ill, we present our results for the quasiparticle DOS with 
and without superconductivity and also break up the re- 
sults into their partial contributions from each of the two 
energy branches separately. We discuss how the pseudo- 
gap alone introduces an asymmetry between positive and 
negative biases in the DOS N{uj) and how superconduc- 
tivity overrides this effect and so restores particle-hole 
symmetry at small energies of order the superconducting 
gap energy A^^. At higher energy, asymmetry remains. 
We also show that modifications in the DOS introduced 
by superconductivity which are confined to a range of a 
few times A^^. in the case of a Fermi liquid extend instead 
over the range of the pseudogap energy scale in YRZ. In 
Sec. IV, we consider decomposing the total quasipar- 
ticle DOS into partial contributions from fixed angle 9 
measured from (tt, tt) in the upper right quadrant of the 
Brillouin zone. These partial distributions can display 
several structures. Nevertheless, upon examination we 
are able to define for each direction 9 an energy gap asso- 
ciated with a specific peak in the partial DOS. This peak 
which is taken mainly, but not always, to correspond to 
the smallest energy closest to a; = 0, is sometimes the 
superconducting gap peak but can also be a pseudogap 
peak. The energies extracted in this way show consid- 
erable anisotropy between positive and negative biases 
which is a fundamental characteristic of the model used 
here. We find evidence in experiment for this anisotropy 
and provide a comparison with STS data. In Sec. V, we 
discuss the limit of low bias where the DOS is linear in oj. 
We recover the FL result with an extra Gutzwiller factor 
reflecting the strong correlations in the system. Section 
VI contains a summary and conclusions. 



II. FORMALISM 



The spectral function for the coherent part of the elec- 
tronic propagator in the model of Yang, Rice and Zhangi 
takes the form 



Q = ± 



gtix)W^[iu'^^rS{uj-E^) + iv^rSiL.+E^)], 



(1) 

where Eg — ^ {E^Y + ^sc(^) the quasiparticle en- 
ergies in the superconducting state with (u^)^ and (f^)^ 
the corresponding Bogoliubov amplitudes. In Eq. ([l}, 
the (a = ±) are the weighting factors of the YRZ 
theory which involve the input pseudogap Apg(fc) and do 
not change in the superconducting state. They depend 
on the electronic band structure energies as well as 
the Umklapp surface energy Specifically, 
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with 6= - (a + Cfc)/2 and Ek = 
the nonsuperconducting state, the 



A2g(fc). 



In 



energies have two 



branches Ej^ = {^k — (^^) /2 ± E^ and in the supercon- 
ducting state, the Bogoliubov weights are given in terms 
of these by 
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In the YRZ paper, the band energies, taken to 
include up to third nearest neighbor hopping, are 
= ~2t{x){coskxa + coskyo) — At' [x) coskxacoskya — 
2t" {x){co&2kxa -\- cos2kya) — fip, where fip is the chem- 
ical potential adjusted to obtain the correct number 
of electrons, through a procedure based on Luttinger's 
theorem. The Umklapp surface energy is where = 
—2t{x) {cos kxa + cos ky a) equals zero. Here a is the two- 
dimensional Cu02 plane lattice parameter. The form of 
the hopping coefficients t{x), t'{x), and t"{x) are fixed 
in the YRZ model^ and will not be changed in this work 
except to note that to enters as a proportionality factor 
in these hoppings and consequently all our results scale 
by to and so this parameter can be varied at will. The 
Gutzwiller factor gt{x) which appears as a simple mul- 
tiplicative factor in Eq. ([l} provides a measure of the 
remaining quasiparticle strength in the coherent part of 
the Green's function. Along with a second Gutzwiller 
factor (?s(x), it also enters the band structure, through 
t{x),t'{x) and t"{x), which provides narrower bands as 
the doping x is reduced towards the Mott insulating 
state at half filling. Specifically, gt{x) = 2a;/(l -f x) and 
=4/(l + x)2. 



For the input superconducting gap, YRZ take 
A" (x) 

A.c(fc) = ^(cosfc.a-cosV) (5) 

and likewise the same form for Apg(fc) with the gap am- 
phtude Apg(x) replacing the superconducting gap am- 
plitude in Eq. ([S])- Both have the simplest d-wave form 
characterized by the lowest harmonic having the required 
symmetry. The amplitudes in Eq. ([5]) and equivalently 
for the pseudogap are 

A°^{x) = 0.14to[l-82.6(a;- 0.2)2] 

A°g(x) = Mo{0.2-x), (7) 

where both are taken to be proportional to Iq and the 
QCP at which the pseudogap gap becomes nonzero is 
Xc = 0.2, where the superconducting gap is also taken to 
have its optimum value. This last condition can easily 
be relaxed to have the maximum gap at 0.16 instead (as 
in experiment). 

In Fig. [TJ we show how the Fermi surface contours 
evolve with doping. At the QCP {x ~ 0.2) there is no 
pseudogap in the YRZ model and the Fermi surface is the 
usual large open surface of Fermi liquid theory (FLT). As 
X is lowered into the underdoped regime, the Fermi sur- 
face reconstructs and Luttinger contours of zero energy 
emerge. They correspond to either = or — 0. 
The solutions for the first case give the hole pockets cen- 
tered on the nodal direction and these exist for 9 — tt/A 
to 9h as indicated in the figure. For the case of a; = 0.19 
but not for the other values of doping shown, there is an 
additional electron pocket (E"^ = 0) located in the region 
between the AFBZ boundary (dashed red line) and the 
Brillouin zone which extends from 9 = 9e to 9 = 0, where 
9 is an angle measured from the origin (tt, tt) in the right 
hand upper quadrant of the Brillouin zone as shown. In 
both cases for hole and electron contours, solutions to the 
equation E^ = when they exist (for a given angle 9) al- 
ways come in pairs. For some angles however, there is no 
solution at all. When two solutions exist, the backside of 
the Luttinger pocket closest to the AFBZ boundary (red 
dashed line) has only a small weight as compared 
with the front part which is orientated towards the cen- 
ter of the Brillouin zone for the hole pocket and oppo- 
sitely for the electron pocket. Both of these have weight 
of order one and this is the piece of the Fermi surface 
(FS) which corresponds to the Fermi liquid (FL) when 
the pseudogap goes to zero. The backsides instead go 
into the AFBZ boundary and have weight exactly equal 
to zero in this same limit. As seen in the figure, when x 
moves towards the Mott insulating state, the Luttinger 
hole pocket becomes increasingly short with 9fi moving 
toward the nodal direction and this is the agency whereby 
the metallicity of the material is increasingly reduced. 
The number of states with significant weighting and zero 
excitation energy is reduced. The approach to half fill- 
ing has a progressively stronger detrimental effect on the 
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FIG. 1: (Color online) The reconstructed Fermi surface con- 
tours for three values of doping: x — 0.10 (brown), 0.16 (yel- 
low), and 0.19 (black), shown in the upper right quadrant 
of the square Brillouin zone. The black contours have both 
hole and electron Luttinger pockets (located around the nodal 
and antinodal directions, respectively). The angles measured 
from (tt, tt) which define the end of these pockets are 6^ and 
9e, respectively. The inset shows the superconducting dome 
and pseudogap line as a function of doping defined via A'^^ 
and Apg in the YRZ model. 



dynamics of the charge carriers due to the increased mag- 
nitude of the pseudogap which has opened on the AFBZ 
boundary. 

An important point to note is that the branch Ej^ cor- 
responds to negative energy except for momenta forming 
the hole pocket, while E'j^ corresponds to positive ener- 
gies except for momenta inside the electron pocket. 



III. RESULTS FOR THE DENSITY OF 
QUASIPARTICLE STATES 

In Fig. [21 we show the quasiparticle DOS correspond- 
ing to each band separately in frames (b) and (c), {oj) 
and iV+(a;), respectively, with the total N{uj) — N~{lu) + 
N^{uj) given in frame (a). Here, N^{uj) is the sum over 
the Brillouin zone of gt{x)W^6{uj - E^). The double- 
dashed dotted black curve for N~{uj) gives results for 
X = 0.2 which has no pseudogap and corresponds to 
the usual Fermi liquid band structure. If to this we add 
the contribution from A'^"'"(a;) in frame (c), we obtain the 
usual FL result shown as the double-dashed dotted black 
curve for the total DOS in (a). It extends over an energy 
range of order several and has a van Hove singularity 
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FIG. 2: (Color online) Density of quasiparticle states N{lj) 
in the pseudogap state as a function of uj in units of to- (a) 
shows the total DOS, and (b) and (c) give the partial results 
from the Ej^ and E'j^ branches separately. The shaded yellow 
region in (b) and (c) identifies the contribution of hole and 
electron pockets, respectively. 



FIG. 3: (Color online) Density of quasiparticle states N(lu) 
when both the superconducting and pseudogap are included, 
as a function of uj in units of to- Similar to Fig. [2] (a) shows 
the total DOS and (b) and (c) give the partial results for ^ 
and E^ g branches, respectively. 



at an energy slightly above — 0.25to- The energy scale 
on which the DOS can vary significantly, however, is set 
by the band width except for the rapid variation around 
the van Hove singularity, but we will not emphasize this 
aspect. Returning to the frames (b) and (c) of Fig. [21 
the areas shaded in yellow for emphasis correspond, re- 
spectively, to the contributions for hole and electron Lut- 
tinger pockets. Note in particular that the shaded region 
for N'^iuj) exists only for dopings near optimum which 
is the only regime for which electron pockets appear. As 
the doping is decreased towards the Mott transition, a 
gap forms in this band above uj = and there is a sharp 
rise in DOS at some finite frequency above which states 
are seen to pile up before iV+(a;) returns to a value closer 
to its no pseudogap value. The energy associated with 
the sharp rise in iV+(w) can be identified as an effective 
pseudogap value, it is different from the input pseudo- 
gap although it is of the same order in magnitude. In 
fact, the net effect of the input pseudogap can be a de- 
pression which can even fall at negative energies as can 
be seen clearly in the composite curve shown in the top 



frame (a) for x = 0.18. The DOS coming from the neg- 
ative branch Ej^ also shows a clear effect of pseudogap 
formation when compared with the double-dashed dotted 
black curve. We see a shift of the main van Hove singu- 
larity (present in the FL case) to lower energies, followed 
by a depression at higher energies beyond a second rela- 
tively smaller van Hove structure. The position in energy 
of this second structure can be taken as a second measure 
of an effective pseudogap. 

In the upper curve for the total DOS, the two scales 
identified as due to the pseudogap combine to give a dip 
in the FL DOS which initially, for small pseudogap, is 
confined to the region of negative energies and is not 
at the Fermi surface as is often assumed in phenomeno- 
logical models of the pseudogap. This dip grows with 
decreasing x both in range over which it extends and in 
depth. For small x it ranges over a large energy region 
and becomes a dominant feature in the DOS which is also 
greatly reduced around the Fermi energy w = 0. It is im- 
portant to emphasize that the pseudogap can introduce 
significant anisotropy between positive and negative val- 
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FIG. 4: (Color online) Comparison of the total DOS N(uj) 
versus w in units of to for doping x = 0.16. Superconductivity 
is included in both curves. The dashed blue curve includes 
the pseudogap. 



ues of uj beyond the relatively mild particle-hole asymme- 
try of the starting FL band structure. This is an essential 
element of the YRZ model which would not arise if the 
pseudogap opened on the Fermi surface rather than on 
the AFBZ boundary. In this regard Fig. [3] is particularly 
relevant. It shows our results when, in addition to a pseu- 
dogap, we include a superconducting gap. Because this 
second gap opens on the Fermi surface it produces a new 
total DOS which is much more particle-hole symmetric 
about ijj — Q than was the underlying pseudogap-only 
DOS, for energies of order of the superconducting gap. 
Before emphasizing this important point further, we note 
that the opening of the superconducting gap has pushed 
some spectral weight in (uj) to lower negative energies 
due to the Bogoliubov coherence factors characteristic of 
Cooper pairing. The peak at the gap at negative energies 
is quite significant in magnitude. 

In Fig. m we compare results of calculations for x — 
0.16 when only a superconducting gap is present (solid 
red curve) and when in addition there is also a pseudogap 
(dashed blue curve). For the solid curve, the supercon- 
ducting coherence peaks fall symmetrically at w = ±0.11 
even though the underlying FL DOS has a van Hove sin- 
gularity which appears only at negative bias. There is 
as well a second peak at w = ±0.23 which is the nor- 
mal state van Hove singularity slightly shifted by the su- 
perconductivity and now appearing at both positive and 
negative biases although its positive bias image is greatly 
suppressed in amplitude. More surprisingly for the pseu- 
dogap case, the symmetry between positive and negative 
w remains in that there are peaks at a; = ±0.08 and 
±0.12 as well as at w = ±0.22 and even at w = ±0.30. 
But of course the magnitude of each peak in a given pair 
can be very different. The asymmetry prominent in the 
nonsuperconducting state has been greatly suppressed by 
the onset of the superconducting gap. Finally note that 




FIG. 5: (Color online) Energy scales extracted from total 
quasiparticle density of states as a function of doping x. Black 
solid and dashed curves are for comparison and are the ampli- 
tudes of the input gap Age (a;) and Apg(a;), respectively. The 
solid red squares are taken from the superconducting coher- 
ence peak position when the pseudogap is set to zero and the 
open red squares are cos(2Sh) times the value of the solid red 
squares. With the pseudogap on, the positive and negative 
biases give different energy scales associated with resolvable 
peaks in N{ui). The green triangles are the superconducting 
gap peaks at larger dopings x > 0.16 and the solid blue dots 
are secondary superconducting gap peaks associated with the 
end of the Luttinger hole pocket. For x < 0.16, the filled 
and open triangles are the pseudogap peaks for positive and 
negative bias, respectively. 



at low bias, the dashed blue curve and solid red curve 
agree very well and the introduction of the pseudogap 
has had no effect on this region of energy. We will return 
to this issue later in Sec. V. 

Having just emphasized the restoration of particle-hole 
symmetry through superconductivity, we explore next 
the asymmetry which nevertheless remains in the quasi- 
particle DOS and which can be traced to the pseudo- 
gap. In Fig. \5\ the solid black curve gives the input 
superconducting gap amplitude Ag(.(x) as a function of 
doping from x = 0.2 (the QCP) to x = 0.1, close to 
the end of the superconducting dome for the parameters 
used by YRZ. The solid red squares were obtained in the 
FL calculations based on the underlying large Fermi sur- 
faces assuming Apg(x) = at all dopings. These points 
are taken from the energy of the superconducting coher- 
ence peaks in these calculations and fall slightly below 
the input values for Ag(,(a:). This is expected since the 
peak in the DOS is representative of the extremal value 
of the superconducting gap on the Fermi surface rather 
than in the Brillouin zone and these are slightly different. 
When both the superconductivity and the pseudogap are 
present, the peak structure in N{u!) becomes more com- 
plex. The lowest energy peaks in the top frame of Fig. [3] 
are associated with the superconducting gap. Starting 
with the most underdoped case first, solid blue curve for 
X = 0.1, we see suppressed coherence peaks as compared 
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to their magnitude in the optimum case x = 0.2 (black 
double dashed-dotted curve). Also, the energy at which 
they occur is considerably less than the value of the input 
gap amplitude. We plot these in Fig. [S] as the solid blue 
circles and see that they trace a dome, but in all cases 
fall considerably below the input gap amplitude curve 
(solid black line). These points correspond to the gap at 
the edge of the Luttinger hole pocket and are therefore 
considerably smaller in energy than Ag(.(a;). In fact, we 
show as open red squares, the product of the solid red 
squares times cos(20^) (see Fig. [T]) and these agree very 
well with the solid blue circles. In all cases, the solid 
and open green triangles for positive and negative bias, 
respectively, are the energies associated with the second 
significant or resolvable peak closest to w = 0. These 
points fall very close to the superconducting coherence 
peak energies of the corresponding FL for a; > 0.16 (the 
triangles are underneath the solid squares for x > 0.17 in 
Fig. [5]) and are the same for positive and negative bias 
w, i.e. we have particle-hole symmetry in this case. We 
also conclude that for this range of doping, these peaks 
are the primary superconducting coherence peaks even 
when the pseudogap is present and further they remain 
largely unmodified from the FL Apg(a;) = case. In 
this sense for these cases, the superconducting gap has 
largely overridden the effect of the pseudogap. This all 
changes for x < 0.16. In those cases, the energy of the 
second peak is quite different for positive and negative 
bias and these pseudogap peaks exhibit a great deal of 
anisotropy. This second energy scale is also seen in Ra- 
man scattering^! where the nodal i?2g geometry probes 
the superconducting gap scale and the Big antinodal ge- 
ometry, the pseudogap. As we have seen here, when op- 
timum doping is approached, the two scales are found to 
merge into a single superconducting gap scale. What is 
different, however, is that the DOS peak structure can 
be exploited to get information on the asymmetric effect 
of the pseudogap between positive and negative energies. 
Our results in Fig. [5] show that this can be considerable 
and that the effect sets in quite abruptly around x — 0.16 
when only the hole Luttinger pocket remains. 

Fig. [5] shows results for the DOS at two dopings x = 
0.14 (left frames) and x = 0.18 (right frames). The top 
frames compare results for normal (red dashed curve) 
and superconducting state (black solid curve). The yel- 
low solid shading and blue hatched regions help one to see 
that our calculations conserve the number of states be- 
tween the normal and superconducting case as they must. 
What is clear from the figure and what we wish to empha- 
size here is that when there is no pseudogap, the states 
lost below the superconducting gap are largely recovered 
just above the coherence peak except for a small shift in 
the van Hove (VH) singularity which introduces a slight 
complication that would not be present in models with 
a constant DOS. When the pseudogap is included there 
is further readjustment of spectral weight introduced by 
the transition to the superconducting state. The energy 
range over which there remains significant changes is now 
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FIG. 6: (Color online) The quasiparticle DOS N{uj) with the 
pseudogap off [(a) and (b)] and on [(c) and (d)] for two doping 
values X = 0.14 (left frames) and 0.18 (right frames). Com- 
parison is made between the normal state (red dashed curve) 
which has no superconductivity and the case with supercon- 
ductivity (black solid curve). The yellow solid shading and the 
blue hatched shading help to see the conservation of states. 



set by the pseudogap scale. Here we are not emphasizing 
the slight complication brought about by the existence of 
a van Hove singularity in the underlying band structure 
chosen in the YRZ model. 



IV. DECOMPOSITION OF DOS IN ANGLES IN 
THE BRILLOUIN ZONE 

We turn next to the decomposition of the total DOS 
into partial contributions coming from different angular 
slices in the Brillouin zone. Such a decomposition has 
been considered by Pushp et ali^ in relation to their 
experimental scanning tunneling spectroscopy (STS) re- 
sults. They effectively decompose the total N{ui) by writ- 
ing 

N{u;)^ -N{u;,e), (8) 

where for N{uj,9), they take the explicit analytic func- 
tional form 

N{ij,9)^Rel , "^^'^ \w{0), (9) 

where F is a smearing parameter to be varied to get a 
best fit of the data along with the value of the gap A{6) 
and weighting factor W{9). In this way, a gap scale A{d) 
can be obtained as a function of 9. Inspired by the above. 
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FIG. 7: (Color online) The partial DOS N(uj, 9) as a function 
of ui in units of to for two values of doping (a) x = 0.12 and 
(b) X = 0.18. From top to bottom, the curves are shown from 
61 = 7r/4 to 6* = in steps of tt/IOO. 



we decompose the full density of states 

7VH = ^A(fe,c.), (10) 

k 

where A{k^Lj) is given in Eq. ([1]), into an integration 
over the magnitude of momentum k for fixed angle 9 
in the Brillouin zone. This gives N{uj,6) directly and 
results are presented in Fig. [71 Frame (a) is for doping 
X = 0.12 and (b) for x = 0.18. Twenty-six values of 
6 are shown between 45° and 0° as measured from the 
(tt, tt) point of the Brillouin zone (see Fig. [IJ. The nodal 
direction (top curves) show small coherence peaks due to 
superconductivity which are centered about a; = and 
are particle-hole symmetric. At higher positive energies, 
we see that there are two more pseudogap peaks while at 
negative bias, a van Hove singularity is seen. 

In Fig.[8l we show results for E'^ and as a function 
of absolute value of momentum |A;| = fc measured from 
(tt, tt) for several values of angle 9 measured similarly as 
the angles shown in Fig. [T] and used for Fig. [T] Doping 
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FIG. 8: (Color online) The variation of the energies E'^ and 
with magnitude of momentum k = |fc|, measured from 
(tt, tt), for various fixed values of the angle 9 defined in Fig. [T] 
as indicated in the figure. The inset gives the energy of the 
pseudogap Dirac point as a function of doping x. Results for 
our complete numerical calculations from curves as shown in 
Fig [7] are compared with the approximate formula -Bdirac = 
{U" ~ Mp)/2. 



was set at a; = 0.14. The pseudogap Dirac point falls at 
9 = tt/A and k = 2.45 in units of 1/a (solid black curves) 
and this corresponds to E'^ = E^ — 0.24 in units of the 
unrenormalized nearest neighbor hopping parameter tg. 
As the angle 9 is increased away from the nodal direction, 
the energies E^ and E^ no longer meet but split and 
there is a clear gap between them. Both move down in 
energy but not by the same amount. Also, the maximum 
and minimum corresponding to a given pair of curves 
do not fall at exactly the same value of k. In the above 
sense, the two van Hove singularities at issue which define 
the pseudogap peaks in N{9,uj) of Fig. [7] are not entirely 
symmetric. The energies of the Dirac point can easily 
be traced as a function of doping. It corresponds to the 
point E'^ = E^ which occurs for S,k + ^1 — 0, i-e. E^ = 

^Jil + A2g(fc) = for 61 = 7r/4, and the details of the 

dispersion curves determine the critical value oi k = kc- 
For X — 0.14, this kc = 2.45 in Fig. [8]which is somewhat 
larger than the value which corresponds to kx = 

ky = tt/2. For this latter point, the dispersion curves 
of YRZ are particularly simple and only third neighbor 
hopping survives. Taking Ef. — but using k^ = ky = 
7r/2 to evaluate the remaining expression, provides an 
approximate estimate for the energy of the Dirac point 
^'dirac as [At" (x) — /^p]/2. This estimate is shown in the 
inset of Fig. [8] to be quite adequate and shows how i?dirac 
moves towards zero energy as x decreases, as is also seen 
in Fig. [7] where the pseudogap Dirac point is seen as the 
closing of the pseudogap coherence peaks at 9 = 7r/4. 



Returning now to this Fig. [3 we note that as is de- 
creased towards the antinodal direction, it remains pos- 
sible to identify unambiguously a lowest energy peak and 
the magnitude of the energy at which these peaks occur 
is recorded on Fig.|9]as a gap in units of for each angle. 
For optimum or even near optimum doping, the resulting 
curve is close to a simple cos(20) curve and represents the 
superconducting gap. However, it should be noted that 
the curves for x = 0.18, 0.17 and 0.16 all show an ad- 
ditional structure between 10° and 15° and the overall 
curve does deviate visibly from a simple cos(26'). This is 
not surprising since near x < 0.2, the Fermi surface, as 
we have shown in Fig. [1] has reconstructed from the large 
open surface of FLT to Luttinger surfaces (holes about 
the nodal region and electrons about the antinodal point) 
and the partial DOS (w, 0) has some knowledge of this 
fact. Consequently, the extracted energy scale is not a 
pure superconducting gap scale. Nevertheless, the dis- 
tortions from a pure cos(20) are not large and the curves 
show almost perfect particle-hole symmetry (compare the 
top and bottom frame of Fig. [HI for x > 0.16). 

The situation is quite different and more interesting 
for low dopings approaching closer to the Mott transi- 
tion. In that case, only a Luttinger hole pocket remains 
as shown in Fig. [T] When superconductivity is not con- 
sidered, there are zero energy excitations along the Lut- 
tinger contours which are real Fermi surfaces and these 
are gapped by superconductivity. But for angles smaller 
than Oh, no true Fermi surface exists and consequently, 
the peaks seen in Niuj^ 0) [Fig. [Tlja)] no longer have their 
origin in Asc(fc) but are related rather to the pseudo- 
gap Apg(fe). This can be easily traced in Fig. (TJa) for 
X = 0.12. In the nodal direction, the peaks nearest w = 
are the superconducting coherence peaks and these be- 
come more prominent as the gap opens up. But eventu- 
ally, particularly on the positive bias side, they start to 
lose intensity while at the same time, the pseudogap peak 
at higher energy remains quite intense and it is this peak 
that must eventually be taken if we are to characterize 
N{ijj, 9) for smaller values of 9 with a single energy scale 
as we are doing here. It is clear that we need to jump 
from one energy scale to the other. On the negative bias 
side, however, the progression with decreasing 9 remains 
smooth. We always take the large intensity peak closest 
to the origin w = 0. These facts are reflected in our Fig.|9l 
In the top frame for x < 0.14, we see two very distinct 
gap scales: a superconducting one for 9 near 45° and a 
pseudogap one for 9 going towards zero. The curves have 
a kink at an angle corresponding to 9 = 9h, the end of 
the Luttinger pocket, but otherwise they show a rather 
smooth behavior. On the other hand, for positive biases 
(lower frame of Fig. ^ there is a clear jump in the curves 
as we transfer from superconducting to pseudogap scale 
and this is followed by a progressive drop to lower values 
as 9 decreases towards zero in sharp contrast to the top 
frame for negative a;. 

To make clearer that our results for the frequency de- 
pendence of N{uj, 9) cannot be fit by the simple for- 
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e (deg) 

FIG. 9: (Color online) The energy gap read off the curves 
for N{uj,9) of the type shown in Fig. [7] as a function of 9 
in the Brillouin zone measured from (tt, tt). The top frame is 
for negative bias (occupied states) and the bottom for positive 
bias (unoccupied states) . Note the isotropy between these two 
sets of data for angles around 45° or 7r/4 (nodal direction) and 
the large anisotropy at small angles (antinodal). 



mula of Eq. ([3]), we have used this formula along with 
the gaps presented in Fig. [9] to recalculate a total DOS 
iV(a;) according to Eq. ([5]) with weights W{9) set equal 
to one for simplicity. Numerical results are presented in 
Fig. [TOlb'). For ease of comparison, we have reproduced 
in Fig. [TUKa) the full DOS already presented in Fig. [3]on 
which our angular decomposition is based. For small val- 
ues of oj, we see a great deal of agreement between these 
two sets of figures. Near optimum doping only the su- 
perconducting gap scale is prominent while for the highly 
underdoped regime both superconducting and pseudogap 
scales are clearly seen. At larger energies important dif- 
ferences arise and these have their origin in the failure 
of Eq. ^ with a single energy scale to capture all the 
details present in the partial densities of states of Fig. [T] 
In particular, the van Hove singularity seen on the nega- 
tive bias side in the top frame is completely missed in the 
lower frame. While this figure speaks to the limitation in 
the method used to extract energy gaps from STS data^S 
it also shows that important qualitative and even semi- 
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FIG. 10: (Color online) Comparison of (a) the original DOS 
N{ij) from Fig. [3] with (b) the quasiparticle DOS Ncb{ijj) ver- 
sus uj in units of to as defined by Eqs. ((8])-(|9ll. Both (a) and (b) 
refer to the state with superconductivity and the pseudogap 
present. 



quantitative information can be extracted in this vi^ay. 
Some details are certainly missed but other important 
features are quite prominently seen such as the pseudo- 
gap scale and its asymmetry. One could improve the 
agreement betvsfeen the top and bottom frame in Fig. [TU] 
by including weights W{9) or multiplying Nog{uj) by an 
envelope function, which partially accounts for the pres- 
ence of the van Hove singularity present in our model 
band structure, but this is not our aim here. 

In Fig. 1111 we provide a direct comparison between our 
results for the asymmetry between positive and negative 
bias of the derived angular dependent gap and the data of 
Pushp et al.— for the same quantity. Plotted is the ratio 
of A_bias/A+bias from our Fig. inifor x = 0.12 (squares) 
and 0.14 (circles). Pushp et al<^ show several curves 
in their figure S3 for a range of data taken from several 
spots on their UD58 sample. As the negative bias data is 
fairly uniform, we took points from along the trend of the 
data. For the positive bias, there is a range of values for 
each angle and so we took the highest and lowest values 
to form separately the ratio with the negative bias data. 
This is shown as the blue triangles and green diamonds, 
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FIG. 11: (Color online) The ratio of A_bias/A+bias versus 
angle, taken from Fig.[9]for x — 0.12 and 0.14. The theoretical 
work is compared with the data of Pushp et al.~ which is 
bounded by the diamonds and triangles (see text for details). 
Inset: a comparison between the theory and experiment for 
the normalized conductance (thick black curve) of Pushp et 
al. and a DOS curve from our calculations (thin red curve). 



respectively. As we wish to concentrate on the asym- 
metry, only the data for angles less than 30° are shown. 
Due to the expected symmetry for angles in the nodal re- 
gion, the data ratio towards the nodal region have been 
normalized to one to facilitate comparison with theory 
even though the ratio in the data was slightly greater 
than one. Our procedure is not entirely rigorous and in 
the hands of the experimentalists there might be a more 
accurate analysis of the result. Nevertheless, while the 
individual curves for A-|_ and A_ are quite different in 
shape as compared with our theoretical results, we find 
that the ratio shows the same qualitative trend as the- 
ory. There is a significant dip around 15-20° followed 
by a rise as the antinodal region is approached. This 
comparison provides evidence that the pseudogap forms 
asymmetrically about the Fermi surface as in our model. 
In the inset to Fig. [TI] we show results for the DOS in 
the resonating valence bond spin liquid compared with 
UD35 data from Pushp et al,^ (their figure S4). We 
used X = 0.11 to match the Tc reduction from optimum 
and have added broadening and a linear incoherent back- 
ground to provide a better fit to the data. This linear 
background does not alter the structures due to the two 
energy gaps in the model. The agreement is good in 
the low energy region shown and it again clearly reveals 
asymmetry between positive and negative biases. 



V. ZERO FREQUENCY LIMIT 

We finally turn to the slope of the DOS at w 
shown in Fig. [T^ a) . In YRZ theory, the superconducting 
Dirac point [shown schematically by Fig. [T^ b)]. on the 
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FIG. 12: (Color online) (a) The slope of N{uj) normalized to 
that for optimal doping versus x compared with the data from 
Pushp et alM'. Inset: N{lli) for a; = 0.14 with both the su- 
perconducting gap and pseudogap present showing the linear 
behavior at low u) in agreement with Eq. (|lip which is shown 
as the red dashed line, (b) A schematic of the superconduct- 
ing Dirac cone as it grows in energy out of the nodal region on 
the heavily weighted side of the Luttinger hole pocket. The 
heavy solid black curve is the Luttinger contour for x — 0.14. 



heavily weighted side of the Luttinger hole pocket in the 
nodal direction, is not affected by pseudogap formation. 
For small energies a;, it is clear from the cone shown in 
Fig. [T21 that the only excited states available are those 
near the bottom tip of the Dirac cone. We can easily 
show that, as for an ordinary FLp^S, the slope is given by 



9t 



TTVAVF 



(11) 



with N{(jj) = s\uj\. Here, vp is the Fermi veloc- 
ity at the Dirac point and va = |VAsc(fc)|fcjj = 
(Agj./'\/2)| sin(fcFx)|, is the superconducting gap velocity 
at this same point. In the inset of Fig. [T^ a). we show 
results of complete numerical calculations of the DOS for 
the case of a; = 0.14 with both pseudogap and supercon- 
ducting gap included (solid black curve) and compare 
with the results of Eq. dTTI) (dashed red line). We see 



good correspondence. 

Formula ([TT|) is the same as would hold in a FL with 
two very important differences: gt and a possible varia- 
tion of VA in this model. In Eq. pTj) . the Gutzwiller co- 
herence factor gt (x) carries the information on how much 
weight remains in the coherent part of the Green's func- 
tion when correlation effects are included, with the rest 
shifted to incoherent processes. This factor is a crucial 
part of our present approach but does not enter FLT. 
An important consequence of this fact is that it reduces 
the rise in the slope in the highly underdoped regime of 
the cuprate phase diagram as compared with a pure FL 
approach. It provides a factor of 2x/{l + x) while the 
gap velocity is proportional to Tc- Pushp et al»^ found 
in STS data that the slope was nearly constant over a 
significant doping range below optimum. Matching their 
reduction from optimum of the Tc of their samples with 
our superconducting dome in order to determine the rela- 
tion to the doping x used in our model, we plot the slope 
of their data normalized to the value at optimum versus 
X. As there was a range of values in their data for a 
particular doping, we have done our best to indicate this 
as a point with a bracketing bar. Comparison with this 
data (solid black circles in Fig. [T2|) shows that, around 
X ~ 0.10 which is the doping in our calculations that 
corresponds to a reduction of Tc by a factor of 3 below 
its value at optimum, the predicted slope represented by 
the dashed red curve starts to increase while experiment 
does not. This could mean that in reality the Gutzwiller 
factor is a more strongly decaying function of x than the 
one we have used. Alternatively, broadening will have the 
tendency to decrease the slope. However, so far we are 
neglecting another important effect associated with the 
YRZ model. In this model the gap to Tc ratio 2A'^^/kBTc, 
which for simplicity we have fixed at a value of 6 in all 
our calculations, is known to vary importantly with x. 
In very recent work, Schachinger and Carbotte'^'^ have 
solved a generalized BCS gap equation with the pseudo- 
gap and Fermi surface reconstruction fully accounted for 
and have found that this ratio changes from its canoni- 
cal value of 4.3 at x = 0.2 (optimum) to ~ 6.5 or even 
higher towards the end of the dome as the Mott insulator 
and antiferromagnetism is approached. Accounting for 
this reduces the variation in slope between optimum and 
highly underdoped by 50%, bringing our calculations 
much closer to experiment as indicated by the red dotted 
curve in Fig. HW a) . The overall agreement between the- 
ory and the data is very good. It is important to stress 
that gt (x) provides an important factor in Eq. (jlip which 
brings theory much closer to experiment than what one 
would find in a FL approach and this also is true as well 
for the variation of gap to critical temperature ratio. 



VI. SUMMARY AND CONCLUSIONS 

We have computed the total quasiparticle DOS in the 
resonating valence bond spin liquid model^ of the under- 
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doped cuprates. When superconductivity is not included, 
the formation of the pseudogap, which provides a mech- 
anism for Fermi surface reconstruction, modifies the un- 
derlying Fermi liquid DOS in an asymmetric way with 
respect to the Fermi energy (w = 0). For small values of 
the pseudogap just below the critical doping associated 
with a QCP, the resulting depression in N{uj) is confined 
to negative energies. As doping is reduced towards the 
Mott insulating state, the depression deepens, covers a 
larger range in energy and spans positive as well as nega- 
tive biases but its effect remains asymmetric. However, if 
superconductivity is also included, particle-hole symme- 
try is restored at lower energies of order A'^^{x), although 
beyond this range the asymmetry associated with the 
pseudogap remains. This effect is traced to the fact that 
the pseudogap is associated with the antiferromagnetic 
Brillouin zone rather than the Fermi surface where the 
superconducting gap opens. 

One can trace peaks in the total quasiparticle DOS 
which are, at optimum doping and just below, coher- 
ence peaks due to superconductivity but these evolve for 
x < 0.16 into pseudogap peaks which are distinctly dif- 
ferent for positive and negative biases. This asymmetry 
is an intrinsic part of the YRZ model and can be used 
to test its validity. We are also able to trace a second 
set of peaks associated with the superconducting gap 
which however originate from the end of the Luttinger 
hole pocket at an angle 9h in the Brillouin zone. In the 
heavily underdoped region of the phase diagram, these 
peaks are the only ones that can be identified as due 
purely to the superconducting gap. These signatures are 
relatively weak, however, as found in, for example, the 
experiments of Boyer et ali^, in comparison with those 
seen in an ordinary d-wave BCS superconductor. As a 
function of doping, the energy corresponding to these co- 
herence structures follow the dome associated with the 
critical temperature although the scale of the effective 
gap involved is reduced. 

When the total DOS is decomposed into partial contri- 
butions from each angle 9 in the Brillouin zone separately, 
N{uj,d), it is found that at some angles N{uj,6) can show 
complex behavior as a function of energy to and even dis- 
play several peaks. Such distributions cannot be modeled 
by the simple function of Eq. Nonetheless, we find 
that we could identify, in a fairly unambiguous fashion, 
a single peak in these distributions which was both close 
to the origin w = and had significant amplitude. This 
provides a single energy scale for each direction 9. At 
angles near the nodal direction, this scale is associated 
with the superconducting gap but as 9 is decreased it 
progresses into a peak associated with the pseudogap en- 



ergy. In some cases the transition from one scale to the 
other is smooth while in others it can be abrupt and also 
somewhat ambiguous. What is found is that the super- 
conducting gap peak at low energies rapidly loses inten- 
sity while the pseudogap peak at higher energies is quite 
intense and so dominates over the lower energy structure 
and thus must be chosen as characteristic of a single gap 
scale for this angle 9 in the Brillouin zone. This gap 
can however differ considerably in size depending on the 
sign of the bias involved. For energy below the super- 
conducting gap scale, there is little asymmetry but this 
changes radically as the antinodal direction is approached 
in the highly underdoped regime where the pseudogap is 
probed. Our findings have implications for the analysis of 
STS data when one wishes to extract directional informa- 
tion from such experiments.^* Comparison with the data 
of Pushp et al. for both N{Ld) and A_bias(^')/A+bias(^') 
shows indeed an asymmetry in the DOS which is in qual- 
itative agreement with the YRZ model. 

Another result is that the zero energy slope of the DOS 
is importantly reduced from its FL value in YRZ theory 
because of the appearance of a Gutzwiller factor gt{x) 
which is a rapidly decreasing function as x decreases and 
represents the reduction, due to correlation effects, of the 
coherent part of the charge carrier Green's function. This 
factor can partially compensate for the rapid increase in 
slope that would occur due to the appearance of the gap 
velocity va in the denominator of the expression for the 
slope in ordinary BCS theory. However, this factor is 
modulated by a rapid increase in the value of gap to crit- 
ical temperature found to be the direct consequence of 
an increase in the pseudogap in the theoretical work of 
Schachinger and Carbotte.'^*' Without these two effects 
there would be a serious conflict between theory and the 
experimental results of Pushp et al.^* who find a slope 
which remains fairly constant in the underdoped regime. 
While a comparison with the existing data and YRZ pre- 
dictions gives qualitative agreement, there are indications 
that the Gutzwiller factor gt{x) and/or the gap to Tc ra- 
tio may in fact vary more rapidly with decreasing x than 
indicated by theory. 
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